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LONG-TERM  GOALS 

One  of  the  major  efforts  in  the  atmospheric  sciences  has  been  to  develop  and  implement  high- 
resolution  forecast  models  and  to  improve  their  parameterization  of  unresolved  physical  processes 
(boundary-layer  transport,  cloud  microphysics...).  For  the  last  three  decades,  the  relatively  pessimistic 
predictions  of  Lorenz  (1969)  about  the  predictability  of  small-scale  (i.e.,  mesoscale)  atmospheric 
features  have  been  largely  ignored  as  routine  weather  forecasts  were  conducted  at  increasingly  fine 
scale.  Recent  research  suggests  there  are  nevertheless,  significant  limitations  to  the  predictability  of 
mesoscale  atmospheric  circulations.  Our  goal  is  to  develop  an  understanding  of  the  predictability  of 
such  circulations  in  forecasts  generated  by  state-of-the-art  high-resolution  mesoscale  models. 

OBJECTIVES 

Specific  questions  addressed  in  our  research  include: 

1 .  How  commonly  does  rapid  growth  of  initial  errors  occur  in  mesoscale  meteorological  settings? 

2.  How  sensitive  are  these  results  to  different  strategies  for  developing  the  initial  ensemble  spread 
using  the  ensemble  Raima n  filter? 

3.  How  can  ensemble  forecasts  be  best  used  to  identify  and  help  predict  these  difficult  events? 

The  answers  to  these  questions  are  of  direct  benefit  to  Navy  forecasters  using  COAMPS  to  produce 
aviation  and  other  forecasts  of  mesoscale  phenomena. 

APPROACH 

The  P.I.  and  graduate  student  Matt  Gingrich,  together  with  Drs.  James  Doyle  and  P.  Alex  Reinecke  of 
NRL.  are  using  the  COAMPS  model  to  conduct  100-member  ensemble  simulations  of  high  impact 
events.  Under  previous  support  we  considered  downslope  windstorms  (Reinecke  and  Durran,  2009). 
which,  it  had  been  argued,  had  high  mesoscale  predictability.  More  recently,  we  have  considered  the 
prediction  of  lowland  snow  in  the  Puget  Sound  lowlands.  Both  of  theses  weather  phenomenon  have 
exhibited  high  sensitivity  to  initial  conditions  in  the  sense  that  the  spread  within  a  large  initial 
ensemble  (either  70  or  100  members)  grew  vary  rapidly  over  time  scales  much  shorter  than  anticipated. 
Part  of  the  motivation  for  this  effort  is  to  help  inform  the  community  of  the  need  to  move  beyond 
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deterministic  mesoscale  forecasts,  which  despite  ail  the  talk  about  ensemble  prediction,  are  still  the 
backbone  of  military,  civilian  and  private  meteorological  forecasts. 

WORK  COMPLETED 

Our  paper  "Large-Scale  Errors  and  Mesoscale  Predictability  in  Pacific  Northwest  Snowstorms’' 

(written  in  collaboration  with  James  Doyle  and  Alex  Reinecke).  was  published  in  the  Journal  of  the 
Atmospheric  Sciences.  A  second  paper.  "Atmospheric  Predictability:  Why  Butterflies  Are  Not  of 
Practical  Importance.”  by  the  P.l.  and  graduate  student  Mark  Gingrich  has  appeared  in  the  Journal  of 
the  Atmospheric  Sciences.  Mark  completed  his  master's  thesis  entitled  Mesoscale  Predictability  and 
Error  Growth  in  Short  Range  Ensemble  Forecasts”  and  also  presented  these  results  at  the  American 
Meteorological  Society  conference  on  Mesoscale  Meteorology  in  August  2013  in  Portland,  Oregon, 

RESULTS 

The  growth  of  mesoscale-forecast  errors  arising  from  uncertainties  in  initial  conditions  was 
investigated  by  examining  100-member  ensemble  forecasts  of  two  powerful  snowstorms  that  struck  the 
East  Coast  of  the  United  States  in  February  and  December  2010.  Our  most  exciting  and  significant 
results  over  the  last  year  connect  the  behavior  of  the  spread  in  ensemble  errors  due  to  initial  condition 
uncertainty  with  error  growth  in  these  two  complex  state-of-the-art  ensemble  forecasts  with  the 
classical  turbulence-model  theory  of  predictability  proposed  by  Lorenz  (1969). 

Many  previous  studies  have  aimed  to  elucidate  the  error  growth  dynamics  controlling  mesoscale 
predictability.  Important  early  theories  were  based  on  an  inverse  error  cascade  in  classical  turbulence 
models,  where  unresolved  small-scale  errors  gradually  propagated  up  to  the  largest  scales  Lorenz 
(1969).  In  contrast,  the  dominant  paradigm  in  operational  mesoscale  meteorology  has  been  one  in 
which  the  mesoscale  is  assumed  to  inherit  predictability  from  the  synoptic  scale  and  thereby  maintain 
forecast  skill  at  much  longer  lead  times  than  those  suggested  by  turbulence  models  (Anthes  1985,  Mass 
2002),  Nevertheless,  the  results  of  this  research,  along  with  several  recent  studies  (Nuss  and  Miller 
2001,  Reinecke  and  Durran  2009.  Durran  el  al.  2013).  suggest  that  mesoscale  circulations  are  in  fact 
extremely  sensitive  to  small  synoptic-scale  errors.  A  variety  of  situations  have  now  been  documented 
in  which  the  degree  of  synoptic-scale  accuracy  required  to  successfully  forecast  mesoscale  weather 
patterns  at  one-  or  two-day  lead  times  would  be  quite  difficult  to  achieve  in  practice. 

Recent  studies  have  suggested  that  forecast  errors  amplify  by  projecting  onto  the  most  rapidly  growing 
physical  structures,  the  scale  of  which  depends  on  the  mode!  resolution  and  the  dynamics  of  the  flow 
being  modeled.  Examples  include  baroclinic  instability  on  the  synoptic  scale  (Tribbia  and  Baumhefner 
2004)  and  convective  instability  on  the  small  scale  (Hohenegger  and  Schar  2007).  Linking  instabilities 
with  upscale  error  growth,  (Tan  et  al.  2004;  Zhang  et  al.  2007)  suggested  a  multistage  process  in  which 
errors  originating  on  the  scale  of  moist  convection  are  responsible  for  stimulating  error  growth  at 
intermediate  scales  that  subsequently  spread  to  scales  large  enough  to  influence  baroclinic  instability. 

The  spectral  structure  of  the  ensemble  spread  in  the  East-Coast  snowstorm  forecasts  was  examined  by 
evaluating  the  ensemble-  and  meridional-averaged  total  and  perturbation  kinetic  energy  spectra  on  the 
highest  resolution  5-km,  convection-permitting  inner  grid.  The  ensembles  clearly  captured  the 
observed  k~5n  total  kinetic  energy  spectrum  at  wavelengths  less  than  approximately  400  km  and  also 
showed  a  transition  toward  a  roughly  k'3  dependence  at  longer  wavelengths.  In  contrast  to  the  small- 
scale  initial  errors  assumed  in  several  idealized  studies  of  atmospheric  predictability,  the  initial 
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perturbation  kinetic  energy  of  our  EnKF-generated  ensembles  was  maximized  at  the  largest  scales. 
This  is  consistent  with  previous  investigations  that  relied  on  data  assimilation  to  either  create  pairs  of 
different  initial  conditions  (Bei  and  Zhang  2007,  Mapcs  et  al.  2008)  or  to  initialize  a  large  ensemble 
(Durran  et  al.  2013),  all  of  which  also  found  that  the  initial  perturbation  kinetic  energy  was  maximized 
at  the  largest  scales.  As  discussed  in  Durran  el  al.  (2013),  this  large-scale  maximum  is  likely  a 
reflection  of  both  small  shifts  in  the  structure  of  the  synoptic-scale  waves  and  the  spectral  signature  of 
isolated,  small-scale  disturbances. 

At  least  as  notable  as  the  initial  structure  of  the  perturbation  kinetic  energy  in  our  ensembles  is  the 
nature  of  the  error  growth.  As  shown  in  Fig.  3,  initial-condition  errors  did  not  simply  propagate 
upscale  according  to  an  inverse  cascade.  Instead,  the  initial  errors  began  growing  immediately  at  all 
scales,  and  the  amplifying  perturbation  kinetic  energy  spectra  formed  a  series  of  self-similar  curves 
over  all  wave  numbers  where  the  error  had  not  yet  saturated.  Following  the  terminology  suggested  by 
Mapes  et  al.  (2008),  the  evolution  of  the  perturbation  kinetic  energy  in  Fig.~17  may  be  described  as 
"up-magnitude"  rather  than  "up-scale". 


Fig.  1:  Ensemble-  and  meridionally-averaged  total  (solid  lines)  and  perturbation  (dashed  lines)  kinetic 
energy  spectra  at  500  hPa  shown  every  six  hours  ( line  colors  given  in  the  legend)  for  the  ensemble 
initialized  (a):  12  UTC  4  February  and  (b):  12  VTC  25  December.  Only  those  wavelengths  greater 
than  the  7 Ax  numerical  dissipation  scale  are  shown 

The  error  growth  in  Fig.  1  can  be  contrasted  with  the  error  growth  in  typical  experiments  with  the 
Lorenz  (1969)  model.  Fig.  2  shows  the  evolution  of  the  error  energy  spectrum  for  the  Lorenz  model  at 
nondimensional  times  in  the  interval  [0,1],  The  initial  error  is  entirely  in  the  smallest  scale  in  Fig.  2a 
and  is  white  noise  in  Fig.  2b.  In  both  cases  the  error  at  the  smallest  scale  is  saturated,  and  this  error 
expands  up-scale  as  progressively  longer  wavelengths  become  saturated.  The  errors  remain  small  in 
those  scales  that  are  not  yet  saturated. 
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Fig.  2:  KE'  spectral  density  as  a  function  of  wavenumber  k  at  non-dimensional  times  t=0  (blue),  0.1 
(green). ...,  1.0  for  initial  error  saturated  at  the  largest  wavenumber  (shortest  wavelength )  and  (a)  zero 
elsewhere,  or  (b)  following  a  white  noise  spectrum. 


hours  (line  colors  given  in  (he  legend):  black  curve  shows  the  saturation  spectrum,  (b)  identical  to  Fig 
l  b.  except  that  the  curves  for  the  total  kinetic  energy  spectral  density  at  each  individual  time  are 
replaced  by  their  average  over  hours  12-36  and  plotted  as  the  heavy  black  line. 

The  evolution  of  the  KE'  spectral  density  in  our  improved  dimensional  version  of  Lorenz's  model  (the 
ssLRS  model)  and  the  COAMPS  ensemble  forecast  from  12  UTC  25  December  2010  are  compared  in 
Fig.  3.  Given  the  extreme  simplicity  of  the  ssLRS  model  (only  24  degrees  of  freedom),  the  agreement 
with  the  COAMPS  ensemble  is  surprisingly  good,  with  relatively  similar  orientations  and  growth  of 
the  KE'  density  spectra  toward  the  saturation  kinetic  energy  spectrum  at  all  times  greater  than  6  hours. 
It  should  be  emphasized  that  the  time  scale  was  not  set  directly,  but  rather  is  determined  from  observed 
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values  for  the  legnth  and  energy  scales.  The  good  agreement  has  provided  empirical  justification  for 
the  use  of  the  ssLRS  model  to  study  the  importance  of  small-scale  initial  errors. 

As  originally  noted,  but  then  largely  overlooked  by  Lorenz,  rapid  downscale  error  propagation  that 
also  occurs  in  systems  with  Ar5/3  kinetic  energy  spectra.  Very  small  initial  errors  in  the  large  scales 
rapidly  propagate  downscale  to  the  shortest  retained  wavelengths.  The  errors  in  the  shortest 
wavelengths  saturate,  and  after  a  brief  period  the  subsequent  upscale  error  growth  is  similar  to  what 
would  have  occurred  if  the  error  was  limited  to  the  smallest  scales  at  the  outset.  Since  the  background 
saturation  kinetic  energy  density  is  much  bigger  at  longer  wavelengths,  very  small  relative  errors  in 
the  large  scales  can  have  the  same  impact  on  predictability  as  saturated  errors  in  the  small  scales. 

, IMPACT/APPLICATIONS 

Forecasting  mesoscate  meteorological  phenomena  is  of  importance  to  many  naval  operations, 
including  those  in  coastal  zones,  those  involving  aviation  in  complex  terrain,  and  those  requiring 
information  about  the  structure  of  the  planetary  boundary  layer.  Understanding  the  degree  of 
confidence  that  can  be  realistically  expected  from  fine-scale  deterministic  weather  forecasts  at  various 
lead  times  will  help  meteorologists  and  other  users  assess  the  importance  of  alternative  approaches, 
such  as  ensemble  forecast  systems.  The  possibility  that  small  initial  errors  in  the  large-scale  analysis 
impose  a  practical  limit  on  mesoscale  predictability  is  a  new  paradigm  that  will  provide  a  further 
impetus  wider  adoption  of  the  ensemble  approach, 

RELATED  PROJECTS 

ONR  Award  Number:  N000141410287  to  the  P.I.  builds  on  these  results  by  continuing  the 
investigation  of  mesoscale  predictability. 


REFERENCES 

Anthes,  R.  A.,  Y.-H.  Kuo,  D.  P.  Baumhefner,  R.  M.  Errico,  and  T.  W.  Bettge,  1985:  Predictability  of 
mesoscale  atmospheric  motions.  Adv.  Geophys.,  28,  159-202. 

Bei,  N.  and  F.  Zhang,  2007:  Impacts  of  initial  condition  errors  on  mesoscale  predictability  of  heavy 
precipitation  along  the  Mei-Yu  front  of  China.  Quart.  J.  Roy.  Meteor.  Soc.,  133,  83-99. 

Durran,  D.R.,  P.A.  Reinecke,  and  J.D.  Doyle,  2013:  Large-Scale  Errors  and  Mesoscale  Predictability 
in  Pacific  Northwest  Snowstorms.  J.  Atmos.  Sci.,  70,  1470-1487 

Hohenegger,  C,  and  C.  Schar,  2007:  Atmospheric  predictability  at  synoptic  versus  cloud-  resolving 
scales.  Bull.  Amer.  Meteor.  Soc.,  88,  1783-1793. 

Mapes,  B.,  S.  Tulich,  T.  Nasuno,  and  M.  Satoh,  2008:  Predictability  aspects  of  global  aqua-  planet 
simulations  with  explicit  convection.  J.  Meteor.  Soc.  Japan.  86A,  175-185. 

Lorenz,  E.  N.,  1969:  The  predictability  of  a  flow  which  possesses  many  scales  of  motion.  Tellus,  21, 
289-307. 


5 


Mass,  C.  F,.  D.  Ovens,  K,.  Westrick,  and  B.  A,  Colie,  2002:  Does  increasing  horizontal  resolution 
produce  more  skillful  forecasts?  Bull.  Amer.  Meteor.  Soc.,  83  (3),  407—430. 

Nuss,  W,  and  D.  Miller.  2001 :  Mesoscale  predictability  under  various  synoptic  regimes.  Nonlinear 
Proc.  in  Geophys.,  8  (6).  429—438. 

Reinecke,  P.  A.  and  D.  R.  Durran.  2009:  Initial-condition  sensitivities  and  the  predictability  of 
downs  lope  w'inds.  J.  Atmos.  Sci..  66  (1 1 ).  3401-3418. 

fan.  Z.-M.,  F.  Zhang.  R.  Rotunno.  and  C.  Snyder,  2004:  Mesoscale  predictability  of  moist  baroclinic 
waves:  Experiments  with  parameterized  convection.  /  Atmos.  Sci.,  61  (14).  1794—1804. 

Tribbia.  J.  and  D,  Baumhefner,  2004:  Scale  interactions  and  atmospheric  predictability:  An  updated 
perspective.  Mon.  Wea.Re y*  132  (3).  703-713. 

Zhang.  F.,  N.  Bei,  R.  Rotpnnd.:  C,  Snyder,  and  C.  C.  Epifanio,  2007:  Mesoscale  predictability  of  moist 
baroclinic  waves:  Com  ectjion-pertnitting  experiments  and  multistage  error  growth  dynamics. / 
Atmos.  Set.,  64  ( 1 0),  3579—359*4. 


PUBLICATIONS 

Durran,  D.R.,  P.A.  Reinecke,  and  J.D,  Doyle,  2013:  Large-Scale  Errors  and  Mesoscale  Predictability 
in  Pacific  Northwest  Snowstorms.  /  Atmos.  Sci.,  70,  1470-1487  [published,  refereed] 

Durran,  D.R.  and  M.A.  Gingrich,  M.A.,  2014:  Atmospheric  Predictability:  Why  Butterflies  Are  Not  of 
Practical  Importance.  /  Atmos.  Sci.,  71,  2476-2488.  [published,  refereed] 

HONORS/AWARDS/PRIZES 

"Atmospheric  Predictability:  Why  Butterflies  Are  Not  of  Practical  Importance.”  by  the  P.l.  and 
graduate  student  Mark  Gingrich  is  currently  ranked  as  number  one  on  the  American  Meteorological 
Society’s  list  of  the  "Top  10  Most  Read  JAS  Articles”  over  the  last  12  months. 

(See  list  at  http://journals.ametsoc.org/Ioi/atsc) 


6 


